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Abstract. I review recent developments in the application of effective field theory to 
nuclear physics. Emphasis is placed on precision two-body calculations and efforts to 
formulate the nuclear shell model in terms of an effective field theory. 



INTRODUCTION 

A question that I have been asked many times is " Why use Effective Field Theory 
in Nuclear Physics" ? The simple and somewhat glib answer to this question is that 
the only other option that one has to using an effective field theory (EFT) is to use 
the " Theory of Everything (TOE), string theory or some derivative thereof. All 
other descriptions must be incomplete at some level and when precise predictions 
are compared with precise measurements, differences will become obvious. It is a 
daunting prospect for us (maybe only me) to use the TOE to compute low-energy 
hadronic processes, and in fact, it is quite silly to even consider such calculations. 
After all, we know that processes in QED can be computed to high precision which 
agree with experimental observations, without knowing anything about physics at 
the Planck scale, M p \. 

The renormalizability of QED assures that ultra-violet divergences, arising from 
our lack of understanding of physics at short-distances, can be explicitly removed 
by a few constants (electric charge and fermion mass), allowing observables to 
be related to each other to arbitrary precision. In contrast, EFT's are non- 
renormalizable but are still predictive when a systematic power counting in small 
expansion parameters can be established. Relations between observables at a given 
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precision will involve a finite number of constants that are not dictated by the sym- 
metries of the EFT alone. For instance, the standard model is a renormalizable field 
theory but processes at energies much less than the scale of electroweak symmetry 
breaking can be described by an non-renormalizable EFT of reduced symmetry, 
SU(3) C Cg> SU(2)l <S> Uy(l) — > SU(3) C ® £/ em (l), where weak interactions are in- 
corporated by higher-dimension four-fermi operators. Even with quantum effects, 
the theory provides a systematic expansion of observables in terms of q 2 /M^ and 
m'j/M^r, where nif is a fermion mass, My/ is the mass of the weak gauge bosons, 
and q is the external momentum. 

If one were interested in calculating the cross section for np — > G?7-radiative neu- 
tron capture by a proton to form a deuteron-directly from QCD, then a lattice 
calculation is the only technique available. The lattice calculation will provide an 
unambiguous cross section in terms of the quark masses and Aqcd, or equivalently 
in terms of other hadronic observables. A cartoon of a contribution to np — > d^ in 
terms of perturbative quarks and gluons is shown in Figure 1. Sources of quarks and 




FIGURE 1. A "cartoon" of a contribution to np — » d'y in QCD. 

gluons that have non-zero overlap with the proton, neutron, deuteron and a source 
for the photon would used to generate the amplitude for np — > dj. It is clear that 
a significant amount of work goes into forming the hadronic states themselves, let 
alone computing the interaction terms. In addition, as the deuteron has such a small 
binding energy and hence is quite extended compared to the nucleon, considerable 
effort will be required to generate the deuteron itself. Unfortunately, at this point 
in time the lattice community is not even close to being able to perform this multi- 
hadron calculation. Indeed, the deuteron itself remains to be generated in lattice 
calculations. 2 If all nuclear lengths scales were of order the chiral symmetry break- 

2 ' During the Effective Field Theory workshop to be held at the Institute for Nuclear Physics at 
the University of Washington during the summer of 2000, efforts will be made to estimate the 
computer resources necessary to determine the deuteron binding energy from lattice QCD [1]. 



ing scale A x then (very) naively lattice computations of nuclear observables would 
not be that much harder than computations in the single-nucleon sector. How- 
ever, there are several low-energy length scales that play important roles in nuclear 
physics. Firstly A x , below which a hadronic description makes sense, and higher- 
dimension operators are induced that describe contributions from scales above A x . 
Secondly, the scale of the repulsive part of the nucleon-nucleon interaction, which 
is conventionally modeled by the exchange of vector mesons (far from mass-shell) 
and numerically is of order A x . Thirdly, the mass of the pion, which is much less 
than A x due to its special status as a pseudo-Goldstone boson. Finally, nuclear 
binding energies which are much smaller than one would naively guess. 

If one is interested in this process at energies much less than A x , or the mass of 
the p-meson, but comparable to the mass of the pion, then it should be sufficient 
to use an EFT with only nucleons, pions and photons as dynamical degrees of 
freedom. All contributions from higher mass scales will be encapsulated in the 
infinite number of higher dimension operators that arise in the momentum and 
chiral expansions. Some diagrams that will contribute to np — ► dj are shown in 
Figure 2. Unlike most EFT's that one encounters, higher dimension operators (dim- 
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FIGURE 2. Contributions to np — ► in a theory with nucleons, pions and photons. Dia- 
gram (a) shows a contribution from 7r-exchange alone, while diagram (b) shows a contribution 
from 7r-exchange and from short-distance interactions. Diagram (c) shows a contribution from a 
local, gauge invariant operator not constrained by nucleon-nucleon scattering data. 

6) involving the nucleon field play a central role [2]- [6]. The fact that the deuteron 
is barely bound, with a binding energy of B = 2.2 MeV, requires a fine-tuning 
between pion exchange and short-distance physics. Naively, one would expect a 
binding energy set by f n , the pion decay constant, much larger than one finds in 



nature. Therefore, the class of diagrams shown in Figure 2(b) is not expected to 
be suppressed compared to those in Figure 2(a). In addition, contributions from 
diagrams shown in Figure 2(c) must be included. These arise from an insertion of 
operators that are gauge invariant by themselves, and are not related in any way 
to operators describing nucleon-nucleon scattering. They arise from short- distance 
physics and have a scale typically set by A x . 

Continuing our descent in energy, if one is interested in this process at energies 
much less than the pion mass, m w , then it should be sufficient to use an EFT with 
only nucleons and photons as dynamical degrees of freedom. All contributions from 
mass scales greater than m n will be encapsulated in the infinite number of higher 
dimension operators that arise in the momentum expansion (chiral symmetry is ex- 
plicitly broken, leaving isospin symmetry as the only relic of the flavor symmetries). 
Some diagrams that will contribute to np — > c/7 are shown in Figure 3. It is im- 




FIGURE 3. Diagrams that contribute to np — > in a theory with nucleons and photons. 
Diagram (a) shows a contribution from interactions between nucleons. Diagram (b) shows a 
contribution from a local, gauge invariant operator that does not contribute (at tree-level) to 
nucleon-nucleon scattering. 

portant to realize that nucleon-nucleon scattering described by this EFT uniquely 
reproduces Effective Range Theory (ERT) [7]. However, for all other observables, 
such as those involving electroweak gauge fields, ERT (e.g. [7,8]) is seen to be an 
uncontrolled approximation to EFT (e.g. [9]). 

In the following sections I attempt to indicate the status of EFT descriptions 
in the various energy regimes. Firstly, I will discuss low-energy |p| <C m n pro- 
cesses involving two and three nucleons, focusing on the recent high precision cal- 
culations that have been performed. Secondly, the issues, results and the present 
roadblocks to successfully describing the intermediate energy regime |p|^ m n are 
presented. Finally, efforts to translate the understanding gained in EFT devel- 
opments to many-nucleon systems are described. Such translation is necessary in 



order to achieve the ultimate goal of having a perturbative theory of nuclei that 
faithfully reproduces QCD. 



NP -> L» 7 AT LOW ENERGIES 

During the past year there has been considerable focus placed on the radiative 
capture process np — > d'y. Firstly, it was pointed out [10] that the uncertainty 
in the cross section for np — > d'y contributed significantly to the uncertainties in 
the predictions of Big-Bang-Nucleosynthesis (BBN) of light element abundances. 
This resulted from the lack of data in the energy region important for BBN from 
either np — > d'y or 'yd — > np. Further, available potential model calculations of 
these processes [11] are undocumented and error estimates are absent. Tools have 
recently been developed that allow for a ~ 1% calculation of the cross section with 
EFT [13,14]. This was facilitated in part by realizing that it is advantageous to 
get not only the location of the deuteron pole correct, but also the normalization 
of the deuteron s-state component. This had long been implemented in the meth- 
ods of [15,16], and implicit in the construction of Weinberg [2,17], but was only 
recently implemented in the dimensionally regulated EFT [12]. Finally, there are 
experimental efforts to measure the small isoscalar £"25 and MI5 amplitudes con- 
tributing to np — > dj using polarized neutrons on polarized protons [18]. At the 
second workshop on Effective Field Theory in Nuclear Physics held at the Institute 
for Nuclear Physics at the University of Washington in 1999, Mannque Rho chal- 
lenged the participants to compute the E2s and MI5 amplitudes with the group 
whose predictions are verified experimentally winning a bottle of exceptional wine. 
With such a wonderful prize at stake, many workshop participants redirected their 
efforts to this project. This is now known as the Rho-challenge. 



np — > for Big-Bang Nucleosynthesis 3 

As existing potential model calculations of np — > d'y are undocumented and error 
estimates unavailable, a 5% uncertainty was assigned to the cross section as input 
into BBN codes [19]. It would be somewhat dismal if, after several decades of 
investigations into nuclear physics, the cross section for this process was uncertain 
at the 5% level. However, EFT calculations have demonstrated that the actual 
uncertainty is much less than 5%. EFT is a well defined method of calculation and 
estimates of the uncertainty in a given calculation can be made by considering the 
magnitude of higher order terms that have been omitted. The expression for the 
cross section for np — > d'y valid at the ~ 3% level (with nonrelativistic kinematics) 
is 

4^( 7 2 +|P| 2 ) 3 r - ~ 2 -I 

a ~ 7 3MMP| M* Ml1 +\ X El\\ ■ (1) 



3 ) I thank Gautam Rupak for allowing me to present his results in this section 



The isovector Ml v and Ely amplitudes are 
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where «i is the isovector nucleoli magnetic moment, a\ = —23.714 ± 0.013 fm is 
the scattering length in the channel, r = 2.73 ±0.03 fm is the effective range in 
the channel, 7 = y/M^B is the deuteron binding momentum, is a number 
derivable from the nucleon-nucleon p-wave amplitudes, Zd = 1/(1 — 7 pa) is the 
residue of the 3 Si nucleon-nucleon amplitude at the deuteron pole, and pd is the 
effective range in the 3 Si channel. There is also a contribution from an operator 
that is not related by gauge invariance to nucleon-nucleon scattering, L np . This 
is the coefficient of a local operator involving four-nucleon fields and a magnetic 
photon. For incident neutrons with speed \v\ = 2200 m/s the cross section for 
capture by protons at rest is measured to be cr^ d = 334.2 ±0.5 mb [20]. The value 
of L np is fixed by requiring that the expressions in eqs. (1) and (2) reproduce cr^f d . 
The amplitudes in eq. (2) have been computed to one higher order by Rupak [14], 
including relativistic effects and the appearance of an El counterterm, providing 
a ~ 1% calculation of np dj. Once L np has been fixed, the photo-dissociation 
cross section 7c? — > np can be determined, and is shown in Figure 4. One finds that 
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FIGURE 4. The photo-dissociation cross section for 7c? — > np. The solid line results from eqs.(l) 
and (2) [13] with L np determined by the cross section for cold np — > dj. The dashed lines denote 
the theoretical uncertainty. Rupak has further reduced this uncertainty to below 1% [14]. 



this relatively simple analytic expression reproduces the data very well, once the 
counterterm L np has been determined at a given energy. An idea of the convergence 
of the effective field theory calculation can be obtained from the numerical results 
of Rupak [14] , 



tr(2 MeV) = 0.0218 (1 + 0.6389 + 0.0135 - 0.0053 - 0.0001 + ...) fm 2 



a(20 keV) = 0.1917 (1 + 0.1076 + 0.0001 + ...) fm 2 



(3) 



which are both seen to converge rapidly. The cross section at various energies 

TABLE 1. a(np — > dj) as a function of the nucleon cen- 
ter-of-mass energy, E. The asterisk denotes an input. 



E (MeV) total EFT a (mb) [14] ENDF a (mb) [11] 



1.264 x 1CT 8 


334.2 W 


332.0 


5.0 x 1CT 4 


1.668(0) 


1.660 


1.0 x 10~ 3 


1.172(0) 


1.193 


5.0 x 10~ 3 


0.4982(0) 


0.496 


1.0 x 10~ 2 


0.3324(0) 


0.324 


5.0 x 10~ 2 


0.1081(0) 


0.108 


0.100 


0.06352(0) 


0.0633 


0.500 


0.0341(1) 


0.0345 


1.00 


0.0349(4) 


0.0342 



computed with EFT by Rupak [14] , along with those from the on-line nuclear data 
center [11] are shown in Table 1. As expected the EFT calculation agrees with the 
numerical values from [11] at the ~ 1% level. 



Isoscalar Ml and E2 Amplitudes in np — > 

As mentioned earlier, the Rho Challenge focused on the Mis and E2s isoscalar 
amplitudes that contribute to np — > d'y, so that predictions can be compared with 
the imminent measurement of polarization observables [18]. Two works were com- 
pleted soon after the challenge was issued, one by Kubodera, Park, Min, and Rho 
[21], and one by Chen, Rupak and myself [22]. 

There are a couple of angular distributions that can be measured in n + p — > d'y, 
but if in addition, the polarization of the 7 can be measured one finds that there is 
a different cross section for production of right-handed versus left-handed circularly 
polarized photons. Defining the asymmetry A 1 (9) to be the ratio of the difference 
to the sum of these cross sections, 

A1(e) =??„ [ (P 7 (M1) + Pj{E2)) cos6 + P 7 (E1) sm 2 6] , (4) 

where the P 7 (nX) are combinations of the M1 V) M1 S) Ely and E2 S amplitudes, 
and r) n is the neutron polarization vector. 

The amplitudes and polarizations are computed in two very different ways. In 
[21] EFT wavefunctions are developed with a coordinate-space cut-off, which are 
then used to determine matrix elements of the various electric and magnetic mul- 
tipole operators. Pions appear as dynamical degrees of freedom and determine the 
long-range part of the nucleon-nucleon interaction. In addition, counterterms are 
included via short-distance interactions (e.g. "delta-shell" and others) so that the 



magnetic and quadrupole moment of the deuteron are recovered. A very differ- 
ent construction is used in [22]. The EFT without pions is used and divergences 
are dimensionally regulated. As in [21], the four-nucleon-one photon counterterms 
are chosen to recover the deuteron magnetic and quadrupole moments, to give 
P 7 (M1) = -7.1 x 10" 4 , P 7 (P2) = -3.5 x 10" 4 and a total of P 7 = -1.06 x 1(T 3 
in the forward direction, approximately 2/3 of the experimentally determined value 
of [23] pexpt _ _(i_5 _|_ 9.3) x io~ 3 . Given the large uncertainty in the calcula- 
tion of the Mis amplitude, and the uncertainty of the measurement, the two are 
not inconsistent. P 7 (M1) = —7.1 x 10~ 4 agrees with the results of Burichenko 
and Kriplovich [24] of P 7 (M1) = —7.0 x 10~ 4 from a Reid soft-core calculation, 
but is somewhat less than their zero-range calculation of P 7 (M1) = —9.2 x 10~ 4 . 
However, given the large uncertainty in the Mis amplitude of [22], both values are 
consistent. P 7 (P2) = -3.5 x 10" 4 calculated in [22] agrees well with that computed 
in [21], and therefore these observables do not distinguish between the two EFT 
methods. 

WEAK INTERACTIONS OF THE DEUTERON 

Weak interaction processes involving the deuteron are central to current research 
efforts in nuclear physics. In addition to the accelerator based programs to elucidate 
the flavor structure of the nucleon, such as the SAMPLE experiments [25] at Bates, 
the interactions between neutrinos and the deuteron form the core of our efforts 
to learn about the neutrino and look beyond the standard model of electroweak 
interactions. Both in production, e.g. pp — > de + u e , and in detection at SNO 
(Sudbury Neutrino Observatory), e.g. v^d — > v^np, charged and neutral current 
weak interaction matrix elements between the deuteron and continuum states are 
required. 

Much effort over the past few decades has been put into calculating the produc- 
tion mechanism pp — > de + u e , both from standard non-relativistic quantum mechan- 
ics [26], and from sophisticated potential model techniques [27]. Recently, EFT has 
been applied to this process by Kong and Ravndal [28] and by Park, Kubodera, Min 
and Rho [29] giving elegant expressions and numerical values of the weak capture 
cross section that are consistent with previous estimates 4 . The cross section de- 
pends somewhat on the value of a four-nucleon-one-weak-gauge-boson interaction, 
with coefficient L^a, as defined in [31]. 

The detection reactions, vd — > npu, v e d — > e~pp and V e d — > e + nn had been 
looked at by two groups, Ying, Haxton and Henley (YHH) [32] and Kubodera 
and Nozawa (KN) [33], using sophisticated potential models. The two sets of 
calculations differ at the 5% level, due to different treatments of meson-exchange 
currents (MEC). Recently, Butler and Chen [31] have determined the break-up 
cross sections with EFT. Only L^a needs to be fixed in order to perform a ~ 1% 
calculation. The YYH and KN numerical results can be recovered with different 



4 ) There has also been recent work in [30] that I have so far failed to comprehend. 
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FIGURE 5. Inelastic v{v)d cross-sections versus incident v(p) energy. The solid curves in the 
left graph are KN results [33] while the dot-dashed curves, which lie on top of the solid curves, are 
NLO in EFT with L x A = 6.3 fm 3 . The solid curves in the right graph are YHH results [32] while 
the dashed curves, which also lie on top of the solid curves, are NLO in EFT with L\^a = 1.0 fm 3 . 
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FIGURE 6. The ratio of charged current to neutral current cross sections in vd scattering versus 
incident v energy at NLO in EFT with Lx a — ~ 20 fm 3 (boxes) and 40 fm 3 (triangles). 



choices of Li )A , as shown in Figure 5 5 , confirming that the difference between the 
two potential-model calculations is short-distance in origin. Therefore, to predict 
the break-up cross section with a precision of better than ~ 5%, one has two options. 
Firstly, compute the /3-decay of tritium, and use this to determine the counterterm 
L\ t A in the EFT, or equivalently the MEC's in the potential models 6 . Secondly, one 
can perform an experiment to measure one of the break-up cross sections to high 
accuracy, and thereby extract Li jA , or the MEC's. Such an experiment is currently 
under consideration [34]. An important input into determining if neutrinos are 
changing flavor as they move out of the sun and to the earth is the ratio of charged 
current to neutral current cross sections. Figure 6 shows that this ratio, unlike the 
individual cross sections, is relatively insensitive to the counterterm Li 7 a- 



a ) I thank Malcolm Butler and Jiunn-Wei Chen for allowing me to reproduce their figures 
6 ) This method of fixing the MEC's has already been implemented for pp — » de + v e [27]. 
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FIGURE 7. Re (<5) versus incident momentum for Nd scattering in the J = | + channel. The 
data are taken from the TUNL pd partial wave analysis [44]. The dashed, solid and dotted lines 
are the LO, NLO and NNLO EFT calculations. 

LOW-ENERGY THREE-BODY PROCESSES 

Significant progress has been made in understanding three-body systems with 
EFT [35,36]. A couple of years ago, Bedaque, Hammer and van Kolck 7 showed very 
clearly that low-energy Nd scattering in the J = | + channel could be described by 
EFT using contact interactions alone. One of the more impressive results was the 
calculation of 03, the scattering length in the quartet s-wave channel. Bedaque and 
van Kolck calculated the first three terms to be a3 = 5.01 + 1.0 + 0.32 + ... fm 

2 

where the ellipses denote higher order contributions, that are estimated to be 
±0.1 fm. The calculated as = 6.32 ± 0.1 fm agrees very well 8 with the experi- 
mental a| xpt = 6.35 ± 0.02 fm [39]. The first term had been computed in 1957 by 

2 

Skornyakov and Ter Martirosian [40] while the second term was computed in 1991 
by Efimov [41]. The third term had not been computed before and was determined 
unambiguously from the nucleon-nucleon scattering amplitude [35] . The results of 
extending this analysis to non-zero energy [35] can be seen in Figure 7 and are 
found to agree well with data. Scattering in higher partial waves has been exam- 
ined by Bedaque, Gabbiani and Griefihammer [37] in the theory with only contact 
interactions between nucleons. A comparison between the EFT calculations 9 , so- 
phisticated potential model calculations and data for one partial wave is shown in 
Figure 8. 



7 ) I thank Paulo Bedaque, Hans- Werner Hammer and Bira van Kolck for allowing me to reproduce 
their figure. 

8 ' Subsequent "second generation" potential-model calculations agree with this result [38] . 
9 ) I thank Paulo Bedaque, Fabrizio Gabbiani and Harald Griefihammer for allowing me to repro- 
duce their figure. 
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FIGURE 8. Re (S) versus incident momentum for Nd scattering in the L — 2 quartet channel. 
The dashed, dot-dashed and solid line is LO, NLO and NNLO. Calculations are shown as crosses 
[42], while the phase shift analysis is shown by open squares and diamonds [43]. 

Computations in the spin-doublet channel required much more development as 
local three-nucleon operators can contribute to the scattering amplitude [36]. The 
scale-dependence of the counterterm (more commonly known as the three-body 
force) is quite different from those that are familiar from perturbative field theory. 
Bedaque, Hammer and van Kolck showed that a single momentum independent 
counterterm could absorb all cut-off dependence from the leading operators re- 
summed by the integral equation that describe three-body systems. The observed 
periodicity as a function of scale indicates that if a given calculation is performed 
with a given value of the cut-off, it is possible that the three-body "force" vanishes, 
while for a different value of the cut-off the three-body "force" may dominate. 
Further, the appearance of only one three-body counterterm required to render 
the scattering amplitude scale independent also naturally explains the Phillips line 
found in potential-models. 

As a last comment on the three-body work, the techniques that have been de- 
veloped in the nuclear physics setting are being applied to atomic systems, most 
notably scattering lengths and recombination rates in Bose condensates [45]. 

ISSUES AT HIGHER ENERGIES 

The situation regarding the applicability of EFT's at higher momentum, near 
or above the pion mass, is much less clear. In Weinberg's scheme [2], pion ex- 
change and local four-nucleon operators contribute to nucleon-nucleon scattering 
at the same order in the expansion parameter. A great calculational and concep- 
tual simplification would arise if the exchange of pions between nucleons could be 
treated in perturbation theory in all partial waves, as suggested by Kaplan, Wise 



and myself (KSW) [6] . Unfortunately, one of the more disappointing results found 
during the past year is that the EFT with perturbative pions [6] appears not to be 
converging [46] (also, see earlier work by Cohen and Hansen [47]). Fleming, Mehen 
and Stewart [46] performed an analytic calculation of the NNLO amplitude for 
nucleon-nucleon scattering in the theory with pions and KSW power-counting [6]. 
They found large non-analytic contributions that appear to destroy the convergence 
of the series. In contrast, several computations were performed with Weinberg's 
power-counting [2] for the nucleon-nucleon potential which appears to give converg- 
ing amplitudes. However, the formal inconsistency of Weinbergs power-counting 
remains. Amplitudes are not renormalization (cut-off) independent at any order 
in Weinberg's expansion, however, the cut-off dependence is found to be numeri- 
cally small when renormalized at a typical strong interaction scale. Therefore, at 
this point in time there is no formally consistent, converging, perturbative EFT to 
describe nuclear interactions for momenta of order or higher than the mass of the 
pion. 

To give you an idea of what has been attempted at these higher energies with 
both Weinberg and KSW power-counting, let me show you the results obtained 
for yd — > yd, deuteron Compton scattering. Three-orders in Weinberg's counting 
have been completed for yd — > yd at higher energies [48] . The angular distribution 
of scattering photons at E 1 = 69 MeV is shown in Figure 9 10 . One can see from 
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FIGURE 9. ^ for yd -» yd for incident photon energy E 1 = 69 MeV with Weinberg's 
power-counting [48]. The dotted, dashed and solid curves are LO, NLO and NNLO. Data is 
from [49]. 

Figure 9 that the expansion appears to be converging nicely to the experimental 
values. Similarly, •yd — > yd has been computed to two orders in KSW power- 
counting [50], the results of which can be seen in Figure 10. Very good agreement 
between data and the parameter free-prediction at NLO is found at E 1 = 49 MeV. 
The agreement is somewhat worse at E 1 = 69 MeV, and does not appear to be 
approaching the data in the same way that the calculation with Weinberg's counting 



10 ) I thank Daniel Phillips, Silas Beane and Bira van Kolck for allowing me to reproduce their 
figure. 
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FIGURE 10. Jft for jd Compton scattering at incident photon energies of E 1 = 49 MeV and 
69 MeV determined in [50]. The dashed curves are LO. The solid (dotted) curves are NLO with 
(without) the graphs that contribute to the polarizability of the nucleon. Data is from [49]. 

appears to. It is clear that higher order calculations must be performed once the 
perturbative pions versus non-perturbative pions issue is understood, and further, 
it is clear that more precise data is required at low-energies. It is worth mentioning 
at this point that neither counting schemes, nor any theoretical calculation that 
exists at present comes close to reproducing the recent data at E 1 = 95 MeV [51]. 
There is much work still to be done in this area. 

ON THE ROAD TO NUCLEI 

In parallel to the efforts that I have described in the two- and three-body sectors, 
Haxton and collaborators [52] have been developing techniques to apply the ideas 
underpinning EFT to the nuclear shell model (efforts are ongoing by others [53] 
but I will not discuss their work in this talk). 

Before I discuss this work I wish to show you the results of a relatively simple 
but demonstrative calculation by Phillips [54] (see also [55]). To show how the 
ideas of EFT can be translated into a potential-model mode of thinking, Phillips 
compared the deuteron quadrupole form factor, presented as T 2 o, computed with 
the Nijmegen93 potential [56] with that generated by an effective potential, Vjg(r) 
(where L denotes the orbital angular momentum state) and local quadrupole mo- 
ment counterterm. The effective potential consists of one-pion exchange at long 
distances and a square- well at short- distances, Vjg(r) = V OPE (r) for r > R, and 
Vgg(r) = Vox for r < R. The values of the V 0i l are chosen to reproduce the 
deuteron binding energy, and low-energy nucleon-nucleon scattering for each choice 
of R. Therefore, the long-distance behavior of the "true" potential and effective 
potential are identical. Further, the tail of the deuteron wave-function produced 
by the "true" potential and effective potential are identical. As one is brutalizing 
the nucleon-nucleon interaction at short-distances while preserving nucleon-nucleon 
scattering, it is expected that predictions for other observables, such as electromag- 
netic form factors, will deviate significantly from nature when probing distance scale 
comparable to or less than R, for reasonable values of R. In addition one expects 
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FIGURE 11. T 2 o computed with the Nijmcgcn93 potential and in the effective theory defined 
in the text for three values of the spatial cut-off R. 

to find that the static moments differ somewhat from nature. It was shown in 
[9,57] that such short-distance modifications can be compensated by the inclusion 
of gauge invariant local operators. In the case of the deuteron quadrupole moment 
it is necessary to introduce a four-nucleon-one-quadrupole-photon operator, that is 
in no way related to the operators determining nucleon-nucleon scattering. These 
operators are induced at the chiral symmetry breaking scale, and must be included 
in any consistent calculation, and in fact, their omission is responsible for the dis- 
crepancy between all sophisticated potential model calculations of the deuteron 
quadrupole moment [58] and its experimental value. Phillips choose a value of this 
operator to reproduce the observed deuteron quadrupole moment for each value of 
R, and then predicted T|q in the effective theory, the results of which are shown 
in Figure ll 11 . It is clear from Figure 11 that by fixing parameters in the effec- 
tive theory to reproduce low-energy observables that, even with a nucleon-nucleon 
potential that has been brutalized at short-distances, one can essentially recover 
the "true" form factor over a quite impressive range of momentum transfers. This 
provides a very clear demonstration that T 20 is determined largely by the tail of 
the deuteron wavefunction, the OPE tail of the nucleon-nucleon interaction and 
the deuteron quadrupole moment. For potential model calculations to accurately 
determine T 2 o, they must first recover the deuteron quadrupole moment, and to do 
so they must include a four-nucleon quadrupole operator [9,57]. 

A very similar exploration is ongoing by Haxton and collaborators [52]. They are 
attempting to construct a model-space-independent shell model, and are presently 
focusing on the deuteron to optimize their techniques. The underpinnings of EFT 
are basis independent, and as such should be able to be implemented in both a 
plane- wave basis (PWB) or a harmonic oscillator basis (HOB). For a continuum 
process, obviously the PWB is preferable, but for a bound state it seems reasonable 
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that a bound state basis, such as the HOB, will be optimal. In an HOB basis it is 
natural to integrate out the levels of the HO step by step to allow for an easier cal- 
culation in a reduced model-space. As each HO level is removed, the hamiltonian 
and coefficients of gauge invariant operators are redefined to preserve observables. 
This is a discrete analog of the renormalization group (RG) implemented in the 
PWB. In Figure 12 12 the magnetic form factor of the deuteron Fmi(q) is shown 
versus momentum transfer, where the n = 50 calculation provides the "true" calcu- 
lation. The other curves correspond to the same calculation, but with the insertion 
of the bare Ml operator in each reduced model space. The solid line, however, is 
not just the n = 50 calculation but also the calculation from ALL reduced model- 
spaces when the renormalized Ml operator is inserted and NOT the bare operator. 
Clearly, a discrete RG can be implemented in an HOB. 

One of the advantages of constructing a discrete RG for the nuclear shell model 
would be to greatly reduce the computer time required to compute matrix elements 
in a nucleus with A ^> 2. Presently, efforts are being made to reduce the shell-model 
space for the deuteron down from n = 140, which reproduces the deuteron binding 
energy perfectly (by construction), down to n = 20 or 30 and faithfully reproduce 
all deuteron observables. Part of the effort is to include the short-range part of 
the nucleon-nucleon interaction by local operators. If successful, this program will 
allow for high precision computations of nuclear properties, with greatly increased 
speed. 
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DISCUSSION 



I have tried to give you an overview of a number of important developments of 
the last year or so. At low-energies, high precision calculations, ~ 1% have been 
performed in the two-body sector, and progress is being made toward calculations 
of similar precision in the three-body sector. 

At somewhat higher energies, a formally consistent and converging EFT describ- 
ing nucleon-nucleon interactions is yet to be uncovered. Weinberg's power-counting 
is formally ill-defined, yet gives numerical results that appear to be converging. 
In contrast, KSW power-counting is formally well-defined, yet appears not to be 
converging! I suspect, as do others, that some sort of union between the two 
power-countings may in fact be the correct one, but this is merely speculation. 

Some of the more interesting developments this year were made in the imple- 
mentation of EFT ideas in nuclear many-body calculations. There are indications 
that the EFT techniques may provide a means to compute properties of nuclei 
presently beyond reach. However, more work is required before any conclusions 
can be drawn. 
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